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ABSTRACT 

We present an analytic calculation of the thermonuclear depletion of the light el- 
ements lithium, beryllium, and boron in fully convective, low-mass stars. Under the 
presumption that the pre-main-sequence star is always fully mixed during contraction, 
we find that the burning of these rare light elements can be computed analytically, 
even when the star is degenerate. Using the effective temperature as a free parameter, 
we constrain the properties of low-mass stars from observational data, independently 
of the uncertainties associated with modeling their atmospheres and convection. Our 
results are in excellent agreement with the detailed calculations of D'Antona & Mazz- 
itelli (1994) and Chabrier, Baraffe, & Plez (1996). Our analytic solution explains the 
dependence of the age at a given level of elemental depletion on the stellar effective 
temperature, nuclear cross sections, and chemical composition. These results are also 
useful as benchmarks to those constructing full stellar models. Most importantly, our 
results allow observers to translate lithium non-detections in young cluster members 
into a mo del- independent minimum age for that cluster. Using this procedure, we 
have found lower limits to the ages of the Pleiades (100 Myr) and Alpha Persei (60 
Myr) clusters. Dating an open cluster using low-mass stars is also independent of tech- 
niques that fit upper main-sequence evolution. Comparison of these methods provides 
crucial information on the amount of convective overshooting (or rotationally induced 
mixing) that occurs during core hydrogen burning in the 5-10 Mq stars typically at 
the main-sequence turnoff for these clusters. We also discuss beryllium depletion in 
pre-main-sequence stars. Recent experimental work on the low energy resonance in 
the 10 B(p,a) 7 Be reaction has greatly enhanced estimates of the destruction rate of 
10 B, making it possible for stars with M ^ 0.1M@ to deplete both 10 B and n B be- 
fore reaching the main sequence. Moreover, there is an interesting range of masses, 
O.O85Af <^ M <^ 0.13Mq, where boron depletion occurs on the main sequence in less 
than a Hubble time, providing a potential "clock" for dating low-mass stars. 
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Subject headings: open clusters and associations — stars: abundances — stars: evo- 
lution — stars: fundamental parameters — stars: low-mass, brown dwarfs — stars: 
pre-main-sequence 



1. Introduction 

Due to low Coulomb barriers and purely strong-interaction reactions, most stars are able to 
destroy deuterium, lithium, beryllium, and boron via proton capture prior to reaching the main 
sequence. For low-mass stars (M <; 0.5 M ; the exact mass cutoff depends on the element, see 
§ light element depletion occurs while the star is fully convective (i.e., before the stellar core 
becomes radiative) and still contracting towards the main sequence. In fully convective stars, 
light elements are depleted throughout the star, and the photospheric abundance is a reliable 
indicator of the interior abundance. In the case of lithium, the depletion at the photosphere 
in these stars is readily observable. Deuterium is depleted on a contraction timescale when the 
central temperature is T c ~ (8-10) x 10 5 K, in which case the energy released from consuming all 
the deuterium (D/H ~ 1.6 x 10 -5 in the local ISM, Piskunov et al. 1997) is comparable to the 



thermal energy content in the star. This slows contraction and nearly leads to a "deuterium main 
sequence" (Grossman &: Graboske 1971, 1973; D'Antona & Mazzitelli 1985) at an age of order 
1 Myr. For some masses, the deuterium burning occurs while the star is still accreting ( |Stahler| 
1988a| ). Once accretion and deuterium burning are complete (after 1-10 Myr, depending on the 



mass, Stahler 1988a) the star undergoes traditional Hayashi contraction. It is during this phase 
that lithium, beryllium, and boron are destroyed in most low-mass stars. 

There has been much previous work on light element depletion in gravitationally contracting 
pre-main-sequence stars (Hayashi & Nakano 1963; Weymann & Moore 1963; Ezer & Cameron 
1963; Bodenheimer 1965, 1966; Nelson, Rappaport, & Chiang 1993; D'Antona & Mazzitelli 1994 
[hereafter DM94]; Chabrier, Baraffe, & Plez 1996 [hereafter CBP96]; Chabrier k Baraffe 1997 
[hereafter CB97]). Motivated by lithium observations along the main sequence in open clusters 



(see Soderblom 1995 and Martin 1998 for recent reviews), many have calculated lithium depletion 
for stars of M ^ 0.5 M & both during and after the pre-main-sequence phase (Vauclair et al. 
1978; D'Antona & Mazzitelli 1984; Proffitt & Michaud 1989; Vandenberg & Poll 1989; Swenson, 



Stringfellow, & Faulkner 1990; Deliyannis, Demarque, & Kawaler 1990; Forestini 1994; Martin 



Claret 1996 ). For stars with M ^ 0.5 M@, most of the lithium burns after convection has halted in 
the stellar core. Accurate predictions of lithium depletion then depend on the temperature at the 
bottom of the retreating convective zone. The location of the convective/radiative boundary and 
the amount of mixing across it depend on opacity, treatment of convection, and rotation; proper 
handling of these effects remains an open question. The range of observational results on lithium 
depletion in young clusters also points to these effects playing an important role for M ^ 0.5 M & 
(Martin 19981 ) . 
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In this paper, we analytically calculate the depletion of lithium, beryllium, and boron for 
those stars that are always fully convective during depletion. The interstellar abundances of these 
elements is low enough that the energy released by their fusion cannot halt (nor appreciably slow) 
stellar contraction. Efficient convection throughout the star allows us to represent it accurately as 
a fully mixed n = 3/2 polytrope, for which analytic treatment of element depletion is tractable 
(Bildsten et al. 1997, hereafter B 2 MU , also see Hayashi & Nakano 1963; Stevenson 1991; Burrows 
& Liebert 1993). Our analytic calculation gives a convenient formula for the time and central 
temperature at which elemental depletion occurs, as well as its dependence on stellar parameters, 
such as the mass and composition of the star. Most importantly, our treatment allows us to survey 
the effects of the presently most uncertain aspect of stellar modeling, the effective temperature T e ^ 
(which sensitively depends on the treatment of the photospheric boundary conditions, opacities, and 
the treatment of convection). These uncertainties have motivated a number of calculations (Pozio 
1991; Magazzu, Martin, & Rebolo 1993; Nelson, Rappaport, & Chiang 1993; DM94; CBP96; CB97), 
each employing different treatments of these uncertain factors. We show that for a given T e g and 
M we can reproduce the results of these prior works to high accuracy. Moreover, our approach 
allows the inferred effective temperature of an observed star to be used directly in analyzing light 
element depletion observations, without first placing it on a theoretical evolutionary track. 



A previous application of our theoretical results to low-mass stars in the Pleiades flB 2 MU|) 



confirmed in a model-independent way the 100 Myr minimum age originally determined from 



lithium observations by Basri et al. (1996). As we discuss in § 5.1.2 , applying our method to 
the lithium non-detections in pre-main-sequence stars in the open cluster Alpha Persei (Zapatero 
Osorio et al. 1996) yields a minimum age of 60 Myr for that cluster, independently of the effective 
temperature scale for low-mass stars. We stress that lithium observations of contracting pre-main- 
sequence stars in young open clusters provide accurate and reliable minimum ages. This is very 
important, because for the age range (10-200 Myr) when low-mass stars (0.09 <^ M/Mq <^ 0.4) 
are actively depleting lithium, stars with masses of 4-15 M Q are leaving the main sequence (using 
Maeder & Meynet's 1989 ty^s(M) relations with overshoot of 1/4 of the pressure scale height at 
the convective/radiative boundary). Confirming or denying that these stars live longer on the main 
sequence due to enhanced mixing (from convective overshoot) is a crucial issue. White dwarfs have 
also been found in many of these young clusters and extrapolating backwards in time from the 
cooling age allows one to estimate the mass of the progenitor star. This extrapolation depends 
critically on the presumed age of the cluster (see Reid 1996 for an up-to-date discussion). 

We begin in § ^ by justifying the treatment of the contracting star as an n = 3/2, fully mixed 
polytrope. We also outline two important limitations to our approach, the formation of a radiative 
core and the onset of substantial hydrogen burning. Our analytic calculation of hydrogen burning 
also yields a convenient fitting formula for the stellar radius and surface gravity as a function of 
T e g along the zero-age main sequence (summarized in Appendix |A]). In § [j| we update the proton 
capture rates for the light elements and discuss the possible effects of a sub-threshold resonance in 
proton captures on beryllium. Section |] presents the equation describing light element depletion, 
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our treatment of the correction due to electron degeneracy at the time of depletion, and our method 
for solving the resulting equations. In § [|, we provide detailed fitting formulae for the depletion of 
different light elements and compare our results to previous numerical work. We emphasize that 
our technique makes robust predictions, even when the instantaneous depletion time is much less 
than the contraction time. This is important, since lithium is detectable even when substantially 
depleted. We also apply our results to lithium observations in open clusters. We conclude in § || 
with a brief summary and a few speculations. 



2. The Internal Structure of Pre Main-Sequence Stars 

Detailed numerical evolutionary calculations of contracting pre-main-sequence (and pre-brown- 
dwarf) stars have been carried out by a number of authors (see Burrows & Liebert 1993 for a review). 
As shown by Nelson, Rappaport, & Joss (1986), very low-mass stars (M <^ 0.03M Q ) remain com- 
pletely and efficiently convective even at ages of about 4.5 Gyr. Higher-mass stars (M ^ O.O7M0) 
stop convecting only once they are degenerate enough (and, likewise, cool enough) for electron 
conduction to carry the necessary flux; this takes at least 5 Gyr ( |Stevenson 1991 ). Of those that 



burn hydrogen, stars of mass M <^ O.3M remain fully convective even after they settle onto the 
main sequence, while stars of mass M <^ 0.5Mq remain fully convective at least until the end of 
Li burning (DM94). This critical mass is lower for beryllium and boron burning. As long as 
we restrict our attention to masses below this critical value (derived in § 2.2), full and efficient 
convection prevails during the period of light element depletion. 

The effective temperature of a fully convective star determines the contraction rate and is found 
by matching the entropy in the deep interior to that below the photosphere (see Stahler 1988b for an 
excellent review of Hayashi contraction). For these low-mass stars with T c q <J 4000 K, the opacities 
are still uncertain, and the treatment of convection near the photosphere is still open to debate. 
The differing input physics results in different effective temperatures for the same stellar mass. For 
example, DM94 found that for a 0.2 Mq star, T e ^ ranges from 3350 K to 3640 K depending on 
whether Alexander et al. (1989) or Kurucz (1991) opacities were used in the stellar envelope. On 
the other hand, for the same star, CB97 get T e s = 3250 K using Alexander & Fergusson (1994) 
opacities, a non-gray atmosphere (Allard &: Hauschildt 1997), and a more consistent treatment of 
convection near the photosphere, as opposed to a gray atmosphere used by DM94. In addition, 
the uncertainty in the effective temperature inferred from the observed spectrum can be as large 
as ±125 K ( Kirkpatrick et al. 1993j ), and discrepancies between effective temperature scales from 



different authors can be much higher, up to 500 K (Martin, Rebolo,& Magazzu 1994; Zapatero 
Osorio et al. 1996|) . 



Despite the variations in the exact value of T e R that stem from differences in the input physics, 
most calculations do agree that T e g remains approximately constant during the fully convective 
contraction phase. Thus, rather than rely on a particular choice of input physics to fix T e f[(M, R), 
we choose to view it as a free parameter that is constant during fully convective contraction, and 
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explore its effect on elemental depletion in a model-independent way. 



2.1. Fully Convective Stars 

In a fully (i.e., core to photosphere; see below for a caveat) and efficiently convecting star 
the entropy per unit mass is independent of position within the star and decreases as the star 
contracts. This is true even if the star is partially degenerate. The entropy of a gas comprised of 
nondegenerate ions of mean molecular weight Mi = pNyi/rii and electrons with weight \x e is 

s = const lnF 1/2 (r?) + -— — tt-ij , (1) 

Mi Me \3Fi/2W J 

where the constant term (which contains the entropy of mixing) depends only on the composition, 77 
is the degeneracy parameter (ratio of the electron chemical potential to fee T), F n {rj) is the nth-order 
Fermi-Dirac function 

F M =r —±* , , <2) 

Jo l+exp(x-7/J 



and F 1 / 2 ('n) x p/p-eT 3 ^ 2 (see, for example, Clayton 1983| ). The entropy is a monotonic function of 



the degeneracy parameter, so for constant \i e and Mi equation (JTJ) is invertible, r/ = r](s). Therefore, 
in an isentropic star the degree of electron degeneracy is constant (Stevenson 1991; Burrows &; 
Liebert 1993; [B^MUp . The equation of state is 



P = P ~^k B T, (3) 

MefT 

where 

J_ = (l + 1 g £Vi^) \ ^ /l 1 (l + 0.1938F 1/2 (r?)) 5 / 3 \ 
Meff Me3Fi/ 2 (?7)y ^ e 1 + 0.12398Fi /2 (r/) /' 

The approximation used in the second equality is from Larson and Demarque (1964) and is accurate 
to within 0.02% for 7] < 30. In essence, for a given pressure and density, the temperature is reduced 
by a degeneracy dependent factor that is constant along an adiabat. This means that the adiabatic 
exponents are Ti = r 2 = T3 = 5/3, and the star is described by an n = 3/2 polytrope regardless 
of the degree of electron degeneracy (Hayashi & Nakano 1963; Stevenson 1991; Burrows & Liebert 
1993; p z MUp . Moreover, Meff is constant throughout the star. Accordingly, we have T = T c 0(£) 



and p = p c 3 ^ 2 (£)> where £ = £ir/R is the radial coordinate, and 0(£) is the Lane-Emden function 
for n = 3/2, which crosses zero at £1 = 3.65. The central density and temperature are 

Pc = 8.44 gem- 3 , (5) 

T c = 



\M Q J V R 



™ * 106 (m) (£) m - w 
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In the non-degenerate limit 1//J, e g — > l//ij + l//u e = and the stellar model is completely 

specified. For an arbitrary degree of degeneracy, however, the stellar radius, R, and // e ff are related 
by 

R 7.73 x lfr 2 /M m \ 1/3 



We will use this relation in § U to calculate the depletion. 

In the above discussion we assumed that efficient convection assures that the star is isentropic 
throughout, and the entire star is well-described by an n = 3/2 polytrope. Electron degeneracy 
alone cannot modify the polytropic structure to affect this relation. However, near the photosphere, 
the entropy cannot be constant for two reasons. First, near the photosphere and above it, radiative 
transport must carry most of the energy flux, and hence the entropy in such regions must increase 
outwards. Moreover, immediately below the photosphere, and especially in the hydrogen ionization 
zone, convection tends to be inefficient and superadiabatic, leading to a decrease of entropy with 
radius; see Figure 3 of Stahler (1988b) for an illustration of this behavior. Therefore, the relation 
between the entropy in the interior and that at the photosphere is not trivial. A significant fraction 
of the effort in numerically simulating pre-main-sequence evolution is aimed at disentangling this 
relationship in order to self-consistently calculate the effective temperatures for these stars. 

The effect of these complications on the stellar radius is extremely small. The radiative skin 
in the stars under consideration is tremendously thin, AR/R < 10~ 4 , AM ~ 10~ 10 Mq ( |Burrows 



Hubbard, & Lunine 1989). Similarly, the ratio of the temperature at the ionization Zone (2~ioniz ^ 
10 4 K) to the central temperature (T c ~ 10 6 K) is of order 10 -2 ; the virial theorem then implies 
that the ionization zone is located in the outer 1% of the stellar radius. Moreover, the ratio of 
the pressure at the ionization zone to the central value is Pi on iz / Pc = (Tiomz/T c ) 5 ^ 2 . Since the 
envelope is geometrically thin, the pressure at the ionization zone is Pi on iz ~ GMAM/AnR, while 
the pressure at the center of an n = 3/2 polytrope is P c = 9.7GM 2 /4ttR 4 , so the mass of neutral gas 
in the star is just AM ~ 10 _4 M. Thus, most of the mass (and volume) of the star is fully ionized, 
as well as fully (i.e., core to photosphere) and efficiently convective, and the stellar structure is 
well-described by an n = 3/2 polytrope until radiative heat transport at the stellar center (§ |2.2p 
or Coulomb corrections (§ pM.l| ) become important. The difference between the true radius of the 
star and that given by equations (|5|) - (|7j) is extremely small, and we are justified in neglecting it. 

Except during deuterium burning, the energy output from burning light elements is insufficient 



to halt the gravitational contraction (however, see {2^ f° r discussion of effects of hydrogen burning). 
Thus, we assume that the luminosity of the star, L, is balanced by energy generation due to 
gravitational contraction alone. For an n = 3/2 polytrope, a virial analysis gives 

r t o , 3 GM 2 dR 

L = 47Ti? 2 a S BT eff 4 = ---^-— ; (8) 

this relation holds even if the star is partially or strongly degenerate. As discussed in § ^, we 
assume that T e ^ is constant during contraction. The time t = then corresponds to the onset of 
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quasistatic contraction, which we approximate with a nonphysical initial state of infinite radius and 
contraction rate. Therefore t differs from the chronological age because of the deuterium burning 
phase, the initial radius on the theoretical stellar birthline ( ptahler 1988a| ), and changes in the 
effective temperature during contraction. However, as long as the effective temperature is correct 
when the depletions occur, t will not differ significantly from the chronological age. Integrating 
equation (||), we obtain the stellar radius as a function of time, 

R 17-1 fM\* t 



R Q Testis 1 / 3 VM G 

where T e3 = 2^/1000 K and t 6 = t/W 6 yr. Since the radius of the star is completely determined 
by Fi^i 7 ])-! we combine the above relation with equation (]?]) to obtain a transcendental equation, 



M F$(,)=4.53xlO-(^){^) , (10) 
for the degeneracy parameter as a function of time. 




2.2. The Onset of Hydrogen Burning and the Formation of a Radiative Core 

As the star approaches the main sequence, p-p reactions generate an increasing amount of 
energy. Eventually, the relation for the contraction rate (eq. ||) ; which assumes that the energy 
generation is only due to gravitational contraction, ceases to hold. This condition defines the lower 
mass limit of our calculations' validity for beryllium and boron depletion during pre-main-sequence 
contraction. We calculate this limit by comparing the nuclear luminosity L nuc (i.e., the total rate 
of energy generation due to p-p reactions in the star) to the energy release rate due to contraction 
(eq. ||). Note that stars depleting lithium either burn it before p-p reactions become an important 
source of energy (however, see § |5.1.1| ), or are not massive enough to reach the main sequence. 

The specific rate of heat production ( Clayton 1983|) for the reactions p + p — > 2 H + e + + v c 



and 2 H + p — > 3 He + 7 is 



l.23x lov^A-V/ir^cxpf-^) <t*s-V. (ID 



where Tg = T/10 6 K, p is in cgs units, X is the hydrogen mass fraction, and / scr is the screening 
correction factor (Salpeter & Van Horn 1969). Further energy release from the PPI chain will 
not occur until the abundance of 3 He has increased to its equilibrium value; the timescale for this 
is greater than 1 Gyr at temperatures ~ 5 x 10 6 K ( |Clayton 1983 ). Since the reaction rate (eq. 



pE ) has such a steep temperature dependence, most burning occurs near the center of the star. 
Therefore, we can approximate the screening factor / scr by its value at the stellar center. Moreover, 



if the star is fully convective, it is described by an n = 3/2 polytrope (see §2.1) regardless of the 
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extent of nuclear burning. Under these conditions, the nuclear luminosity is 

Lnuc = J e pp dm « 4.72 x 1(T 2 T C 6 36 e pp (p c , T C )M, (12) 

where the prefactor in the above equation comes from integrating the reaction rate over an n = 3/2 
polytrope and constructing a simple analytic fit to the result. Using equation (|12"D, we have derived 
an analytic relation between the effective temperature along the zero-age main sequence and stellar 
radius and surface gravity (see Appendix [A]) . 

When the nuclear luminosity is much less than the luminosity due to gravitational contraction 
L c (eq. ||), the assumptions leading to the contraction law (eq. ||) hold, and our depletion cal- 
culations are reliable. We indicate in Figures |l] and ||| the central temperatures and ages at which 
L nuc ~ 0.05L C and L nuc ~ 0.50L C (heavy solid lines). For comparison, we also show (triangles) a 
1 Gyr isochrone (roughly, the main sequence) from CB97. Our results for light element depletion 
are only valid in the regions of the plots below the aforementioned lines. As one can see, lithium 
depletion is unaffected by the hydrogen burning, while the beryllium results are only valid for 
M ^ O.IMq (boron depletion will be discussed in § |5.3j). 



Although p-p reactions can slow the contraction rate, they do not alter the n = 3/2 polytropic 
stellar structure. However, for stars of mass M £ O.35M flChabrier fc Baraffe 19971 ) a radiative 



core forms during contraction. Formation of a radiative core breaks the polytropic structure, 
prevents mixing into the depletion region, and thus defines the upper mass limit of applicability of 
our depletion calculations. We estimate the central temperature at which convection ceases in the 
core as follows. 

The gas at the stellar center is convectively unstable whenever the temperature gradient 
d In T/d In P exceeds the adiabatic value 1 — I/T2 (see, for example, Hansen & Kawaler 1994). 
In other words, the maximum rate of energy generation at the center of the star that can be 
balanced by radiative heat transport alone is 



Crad 



dL(r) 



dm 



64vr Go-sb T c 4 / 1 



p c V r 2 



1 - =r , (13) 



where k is the opacity. The minimum rate of energy production e ra d for which the stellar center is 
fully convective is just equation fll3| ) with T2 = 5/3. In pre-main sequence stars, the bulk of the 
energy is released by gravitational contraction, with the specific energy production rate 

P d f P \ _3k B N A TdlnR 

£grav " - P (r 3 -i)dt n {7^) ~ "2^ dT ; (14) 

the total luminosity L (equation ||) is just the above rate integrated over the entire star. As the 
star approaches the main sequence, there is an additional small contribution e pp due to hydrogen 
burning (eq. [11]). We equate the total energy generation rate at the stellar center e = e grav + e pp 



with e rac j and solve for T c . We compute the opacity at the center from the OPAL tables ( [Rogers 



fc Iglesias 1992) and use effective temperatures as a function of mass from CBP96 to compute 
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the contraction rate. The resulting radiative/convective boundary curves are plotted in Figures || 
and |2| (heavy dot-dashed lines). The stars develop radiative cores in the regions of the plots above 
the e = e ra d curves. For 7 Li, the mass, central temperature, and density at which half-depletion 
is coincident with the formation of a radiative core are M1/2 « O.5M , T c w 4.1 x 10 6 K, and 
p c ~ 4.8 gem -3 . This is in excellent agreement with CB97, who find that a radiative core forms in 
a O.5M star at an age of 10 7 yr, when the lithium abundance is 0.52. The corresponding stellar 
parameters for which the formation of a radiative core is coincident with 7 Li depletion by a factor 
of 100 are Mi/ioo ~ 0.44M3, T c ~ 4.3 x 10 6 K, and p c ~ 7.7gcm~ 3 . For 7 Li, our fully convective 
depletion calculations are therefore valid for M <^ 0.45 — 0.5M Q , depending on the desired depletion 
level. 

We should point out that the above results are quite sensitive to the value of the opacity at 
the stellar center. If we write the opacity as K p Kp T KT (with p and T in cgs units), we find that, 
at the central temperature where 7 Li is depleted by a factor of two, k q = 5.66 x 10 21 cm 2 g -1 , 
Hp = 0.347, and kt = —3.166. The mass of the star that becomes radiative upon depleting half 
of its 7 Li content scales with k q and T e $ as M1/2 oc k^ A8 TJ^ g . The sensitivity of observed lithium 
depletions to the interior opacity has previously been considered for higher mass main-sequence 
stars by [Swenson et al. (1990) and Swenson et al. (1994) , who adjusted the interior opacities by 



changing the metal content. Although they were considering stars of approximately solar mass, 
Swenson et al. (1990) found that an increase of 37% in the interior opacity increased the mass of 
a Hyades star at which the lithium was half-depleted by ~ 20%, i.e., M1/2 oc k^ 5 , similar to our 
estimated dependence. 



3. Proton Capture Rates for Light Elements 

We use thermonuclear reaction rates in the form 

N A (av) = Sf scv T 6 -i exp (-^73) cm 3 s" 1 g" 1 , (15) 

to approximate the rates of Caughlan & Fowler (1988, hereafter CF88). For non-resonant reactions 
j = 2/3, while for reactions affected by resonances the value of j has to be adjusted (see below). 
Here, S does not represent the astrophysical S-factor, but is rather, like a, a dimensionless parameter 
in the fit to the reaction rate. Table |TJ shows the values of S and a found by fitting to the rates 
over the range of temperatures (Tq < 6) appropriate for this work. Raimann (1993) recently 
incorporated new experimental results at low energies (~ 11-13 keV) to update the rates for 7 Li, 
adjusting the value of S from 6.4 x 10 10 (CF88) to 7.2 x 10 10 . 

There are two stable boron isotopes with abundance ratios in the local ISM of 11 B/ 10 B= 3.4^ 'g 
( Fcdcrman et al. 1996| ). The 10 B(p, a) 7 Be reaction is resonant, as an excited state of the n C nucleus 
exists just 10 keV above the center-of-mass energy of the p+ 10 B system. As a non-resonant reaction, 



the center-of-mass energy would be at E c = 3.46T6 2 / 3 keV, and the reaction occurs (as we show 



- 10 - 



later) at Tg ~ 6, so that E c = 11.4 keV, and thus this resonance affects the cross section. Rauscher 
& Raimann (1996) recently corrected the CF88 rates for this broad resonance, which enhances the 
rate by at least a factor of 200 for the temperatures of interest. For T 6 < 6, a reasonable fit to 
the reaction rates that includes the resonance is to multiply the CF88 rates by 300{T G /3) 1 '' 2 . The 
entries for this reaction in Table [l] include the resonant correction. 

There are potential problems with the 9 Be (p, a) 6 Li rate as well, since an excited state of 
10 B (E x = 6.56 MeV) lies just 25.7 keV below the energetic threshold. Current knowledge (see 
Ajzenberg-Selove 1988] ) holds that the angular momentum and parity of this state is either J 71 = 2~ 
or J w = 4~, with 4" thought to be the better guess. The angular momentum and parity of 9 Be 
is (3/2)~, while that of a proton is (l/2) + . Hence, if this resonant state has J n = 4 _ , the proton 
must overcome a centrifugal barrier of angular momentum i = 2, and the resonant state will not 
contribute significantly to the reaction rate. But if instead J 71 = 2~ , the resonance can proceed via 
s-wave (£ = 0) and the rate will be greatly enhanced over the CF88 rate (Brown 1997). 

To estimate the effect such a resonance would have, we parameterize the entrance channel by 
the dimensionless reduced width 9 2 ( |Clayton 1983 ). Crudely speaking, 9 2 measures the degree to 
which the compound p + 9 Be nucleus resembles the excited level of 10 B. The cross section is then 
calculated using the Breit- Wigner formula (see Brown 1997] for details) with 9 2 as a free parameter. 
A good fit to the thermally-averaged rate for J n = 2~ (accurate to within 1% for 2 < Tq < 5) is 

N A (av) « 3.53 x hWtT 1 - 367 exp ("^f) cmVV 1 ; (16) 

the total reaction rate is the sum of this equation and the standard non-resonant rate. For the 
temperatures of interest in low-mass stars, Tg <J 6, an s-wave resonance will dramatically boost the 
rates if 8 2 p> 0.01. In this paper, in addition to using CF88 rates for 9 Be, we use equation ( |l6| ) to 
explore the effect of the resonance on beryllium depletion in pre-main-sequence stars. 



4. The Depletion Equation 

The central temperature needed for depletion can be roughly estimated as follows. Light 
elements are depleted when the nuclear destruction time at the center of the star, 

where n p is the proton density, becomes comparable to the contraction time (i.e., the timescale for 
the change in central temperature if degeneracy is mild), 

R /3000i^T\ 4 /O-1M \ /0.6\ 3 / T c \ 3 , , 

«- s -557* = 115 hr-) («) i^w) Myr ' (18) 

The above equation is valid regardless of the degree of degeneracy so long as the star is well described 
by an n = 3/2 polytrope. However, for rough estimates, one can put /j, e g = fj, = (l/fa + l//i e ) _1 



-11- 



in equations fll7| ) and (18). Evaluating idest = ^cont at the stellar center gives an approximate 



condition for the central temperature at the time of depletion, 

-4tz = 32.9 + ln(5/ scr )-31nf^— ) -4ln(-^-) -6ln(^) + ( 6 - j) lnT c6 . (19) 
T V3 ^ v JSCT) V0.1M o y V3000K; \0.6/ 

c6 

The destruction time is a very steep function of the central temperature: once T c reaches a certain 
value, the entire abundance of a light element is depleted on a timescale much shorter than t con t. For 
example, a 0.1M@ star with T e ff = 3000 K burns Li when the central temperature T c ~ 3 x 10 6 K. 
Under these conditions, the reaction rate is extremely temperature sensitive (Na(ctv) oc T 20 ) and 
the central temperature necessary for elemental depletion depends only weakly on mass. The 
sharp dependence of the reaction rate on the central temperature means that our calculations are 
extremely insensitive to small uncertainties in the constitutive physics, such as screening ( |B Z MU ). 



For example, the central temperature at depletion is proportional to (S'/scr) 1 ^ 20 and the time of 
depletion t oc T 3 oc (S/ scr ) 3 / 20 . An uncertainty in the reaction rate or the screening correction by a 
factor of 2 can change the depletion time by no more than 10%. The exact equation for elemental 
depletion (which we derive below) is just a fancy way of writing t con t = idest that includes all 
the prefactors that arise from integrating the destruction rate over the entire star, and from the 
changing degeneracy (i.e., the decrease of /x e g) during contraction. 

At a given time t, consider a mass shell of coordinates (m, m + dm) that contains f(m, t) dm 
grams of a given light element isotope. There are two contributions to the rate of change of f(m, t): 
(1) proton-capture reactions, and (2) convective transport from isotope-rich regions to isotope- 
poor regions. We ignore diffusion, as it occurs over timescales too long to be of any relevance. 
Symbolically, 

df\ fXp(m,t) . ., . OF 



dt J m mn ' dm ' 

where F is the net mass of the light isotope transported per second by convection across a surface 
r(m) at a time t, where mu is the hydrogen mass. The term dF/dm vanishes upon integration 
over the entire star. Defining the mass-averaged value of / by the equation Mf = f dm, we 
integrate equation (|(]) over the star to obtain 



M% = - ( M —p{av) dm. (21) 
dt Jo mu 

The timescale for a fluid element to traverse the radius is t m j x ~ £/v, where I is the total length 
traveled. For a random walk process (Bodenheimer 1965), £ ~ R 2 /Xp, where Xp is the pressure 



scale height, which is of the order of the mixing length. The velocity may be estimated for the case 
of efficient convection, giving t m i x ~ 10-100 yr in the core. Inserting a value of t appropriate for 
the photosphere gives t m j x ~ 10 3 yr, which is still much shorter than the evolutionary time scale of 
10 6 -10 8 yr. Therefore, the star is well mixed and f(m,t) can be replaced by f(t) in the integral. 
Then, upon changing variables to spatial coordinates in the RHS of equation (pl|), we have 



^ In/ = -i^- ( R p 2 {av)r 2 dr. (22) 
dt J m H M 7 
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For simplicity of notation, we drop the bar over /. 

Using the non-resonant form of the reaction rate (eq. []15| with j = 2/3) and the definitions of 
central density and temperature (eqs. || and H), we write the differential equation describing the 
change of mass fraction / as 

where a = aT c g _1 / 3 . The burning rate is very temperature-sensitive because the center of the star 
is much colder than the Gamow energy. This sensitivity restricts the burning to the central region 
of the star, where the temperature and density are quadratic in radius. Specifically, the peak of the 
rate integrand lies at £ pca k ~ 0.042 (55/a) 2 + 0.0067 (55/a). For this same reason we approximate 
the screening factor by its value at the center of the star. We then find an analytic approximation 
to the above equation that is accurate to within 1% for a > 45 by using the first two terms of the 
expansion 9(£) = 1 - £ 2 / 6 + • • • and taking the upper limit to infinity. This yields: 

u f = -o.i8 (*\ fny r^v«_v"/» d - ^ 



^'-^HotAsJ {-W) Sf^ a -^(l-- )e -^. (24) 
A similar development for resonant rates (with j kept as a parameter) yields 

- - (£) (S) 3 ( t) 2 (> - — 1 <*> 

In order to integrate this equation, we need to know the dependence of the central temperature 
parameter a on time, which is found by using the polytrope equation (^) and the degeneracy 
relation (eq. @): 



2 = o.5i3 (mv /3 f^y /s ( a) i/s = 6 . 73 ( (^m Ila 

a \fi cS J \MJ \R Q J \MqJ \ t 6 T* 3 I 



The value of i ? i/2( r ?) is obtained by solving the transcendental equation (|l0|). Then we can use 



equation ( |26|) to integrate the depletion equation (eq. |24| or eq. [25]) numerically from t = and 
In / = to find the element abundance as a function of time. Results of the calculations based on 
this formalism are presented in § |5[ 

In the region of the parameter space where degeneracy effects are not dominant (i.e., to the 
right of the constant degeneracy line (heavy dotted line) in Figure |l|), the above derivation can 
be considerably simplified, and the depletion equation (either eq. |24| or pq| ) can be integrated 
analytically (see B^MU for a simple version of this development). During contraction, the central 



temperature initially increases because the radius decreases (from the virial theorem for a nonde- 
generate ideal gas T c oc M/R). For a nondegenerate gas, l/// e ff = 1 //Ui + 1 //U e , but when degeneracy 
becomes important /i e g decreases, leading to a decrease in central temperature. Explicitly, we have 

d j din f da dlnf / da dR ^ da dn c g\ 
dt da dt da \dR dt d/j c s dt J 
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When degeneracy is not important (roughly, for M ;> 0.2Mq), /i e g- is negligible compared to R. 
Thus, we can rewrite the depletion equation (p5|) in a simple form 

4 



d 
da 



In/ 



6.23 x 10 



14 



x5/ scr a 18 - 3 ^-( 41 /2-3.) 



X \ /0.6\ D /1000KN 4 /M 
27 - 9j 



1 



2a 



This equation may be integrated directly from the initial condition ao 

6 / 1 nnn xc \ 4 / ^ \ 3 



VF 



, /o 

ln 7 



6.23 x io- (*\ /oey/ioooK 



M 



oo to a: 

18 ~ 3 M«), 



where 



9iaJ 



a -(41/2-3i) e -a 



68 - 15j 



r(-(41/2-3j),a) 



(28) 



(29) 



(30) 



and T(x,a) is the incomplete gamma function. For a given W, equation (29) can be solved for 
a(W), which is essentially the central temperature required for the level of depletion W. This is 
the method used in B^MU to compute 7 Li depletion. We find that the two approaches (neglecting 
or including /i e ff) agree very well whenever the approximate formulation, equation (^8|), is valid. 



We note that the treatment of degeneracy and screening in equation fl29|) is approximate in 
the following sense. In addition to neglecting <9a/<9/x e fj when changing variables from t to a, we 
neglected the dependence of fi e g and / scr on a when integrating equation (28). This is justified 
because, when degeneracy is mild, fi e g and / scr vary very slowly with a, while the radius R is a 
strong function of a. Furthermore, the contraction rate of the star is dR/dt oc i? 4 , so the star spends 
most of the time at smaller radii, where the contraction rate is small. Thus, to approximate electron 
degeneracy and screening corrections we just use the final values of /j, e ff(a,M) and f scr (a,M) in 
equation ( p9|) and solve the resulting transcendental equation. We find that the errors introduced 
by this approximation, as compared to integrating equations 



or (25) directly are never more 
than a few percent. For rough estimates that neglect effects of degeneracy and screening, one can 
just set l/Heft = 1/fM + l/fie = l/n and / scr = 1. 



5. Discussion of the Results for Each Light Element 

We have presumed that the contracting pre-main-sequence star is well described as a fully con- 
vective n = 3/2 polytrope of arbitrary degeneracy with no substantial internal energy generation. 
This limits our light element depletion calculations in three ways: 



1. Onset of Hydrogen Burning. Hydrogen burning begins (see § 2^2) before the star settles 
onto the zero-age main sequence and thereby increases the time needed for the star to reach 
a given radius, as compared to equation @. Since we do not include this effect in our 
calculations, we limit our discussion to those stars generating less than 5% of their luminosity 
by hydrogen burning. 
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Development of a Radiative Core. Depletion at the stellar surface is directly related to 
the central temperature only when the star is fully convective. Further depletion will occur 
in the interior parts of the star, but the surface abundance of a light element is then sensitive 
to the temperature at the radiative/convective boundary and not to the central temperature. 



We stop calculating depletion when the core of the star becomes radiative (see § 2.2) 



3. Coulomb Effects on the Equation of State. Depending on the particular element consid- 
ered, Coulomb corrections to the equation of state can be important at the time of depletion. 
These act to break the polytropic nature of the star and increase the central temperature 
relative to the polytropic value (see § |5.1.1|) , 

Figure p] displays our results for central temperatures when the 7 Li abundance varies from a 
half to a hundredth of its initial value (shaded region). We also plot our central temperatures at 
one-half depletion of 6 Li (light dotted line) and 9 Be (light dashed line). (The boron isotopes are 
discussed in § |5.3| .) The same curves, expressed in terms of mass and age, are shown in Figure || 
and motivate the following discussions for each element. 



5.1. Lithium Depletion 

The substantial amount of lithium]^ observations already makes the case that lithium abun- 



dances are a reliable indicator of stellar ages and masses (Rebolo, Martin, & Magazzu 1992; Martin 



et al. 1994| ; |Basri et al. l99§ ; [Rebolo et al. 19961 ; Pppenh eimer et aL 1997| ). In § |5.1.2 we present 



an extended discussion of the observations of the cluster Alpha Persei and derive a minimum age 
of 60Myr. Since this same technique can be applied to clusters in the age range of 10-200 Myr, we 
extensively discuss the application of our results to the lithium observations. 

Previously (in B 2 MU ), we discussed the lithium depletion results and provided convenient 
fitting formula. In addition to providing slightly improved formulae here, we also discuss in detail 
the limitations of our work at the very low-mass end. The intersection of the lithium depletion 
with the onset of a radiative core (see Figures |l| and |2[) confines our work to M ^ O.46M (for this 
mass, the ending of convection in the core is concurrent with depletion of Li to 1% of its original 
abundance). Our method has the advantage of leaving T e R independent of stellar mass, but the 
disadvantage of not including Coulomb effects, which limits us to masses M ;> 0.072 Mq, or, for 
typical values of T c ^ , to stars younger than ~ 200 — 300 Myr. 



2 Throughout this section references to lithium mean 7 Li, as 6 Li always depletes first and is less abundant in the 
local ISM (Lemoine, Ferlet, & Vidal-Madjar 1995). 
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5.1.1. Comparison to Previous Theoretical Work and Limitations at Low Masses 
Based on the results presented in B Z MU, and using a fiducial effective temperature relation, 



T e g oc M 1 / 7 , we obtain the central temperature at the time of one-half lithium depletion, T c ~ 
3.1 x 10 6 K(M/O.1M0) 1//6 . From the polytropic relation (eq. |5)), the central density at this time 
is p c rs 62 gem" " 3 (O.1M /M) 3 / 2 . From T c and Pc we assess how important Coulomb corrections 
to the equation of state will be at the center of the star when the burning occurs. The plasma 
parameter at the time of one-half depletion is 



Accurate evaluations of the equation of state in this regime are known (see Saumon et al. 1996). 
However, we will roughly estimate the size of the Coulomb correction using the Coulomb pressure 
for T < 1 (Clayton 1983), P c = -e 3 /3(7r//c B r) 1 / 2 (1.8p/m J /) 3 / 2 , which then gives 

Pc _, oa 2 /O.lM 



p yxl ° iirj- (32) 

at the time of one-half depletion. Since the Coulomb pressure is negative, the central temperature 
is higher than its polytropic value. 

Our neglect of the Coulomb correction places a lower limit on the mass of a star for which we 
can reliably compute the light element depletion. It is easy to show that, for normal composition 
{X = 0.7), the central temperature as a function of degeneracy is 



Tc = 5 - 63 x 10 { o^tom; ) ^tfto) K - ^ 



'0 

As the star contracts, the central temperature rises until the degree of degeneracy reaches a critical 
value (-^1/2(77) ~ 4 — 6), after which the central temperature begins to decline. The maximum value 
of T c is 

„ / M \ 4 / 3 

T c max = 3.43 x 10 6 K. (34) 

' VO.O75M / v ; 

We note that equation fl3"4| ) fully incorporates the effects of partial electron degeneracy, while a 
similar relation for T Cimax using an approximate treatment of degeneracy is derived in Stevenson 
(1991) and Burrows & Liebert (1993). The maximum temperature as a function of mass (for 
M < O.IMq) is denoted with a thick dashed line in Figure |l|. The intersection of this curve 
with the lithium depletion line defines the lower mass limit of validity of our lithium depletion 
calculations. 

For high-mass stars that deplete lithium long before the central temperature reaches its peak, 
our neglect of Coulomb corrections slightly delays the time of depletion. This is because the 
Coulomb corrections allow the star to reach a given T c at an earlier time and thereby deplete 
earlier. Roughly, since the time to reach a given central temperature scales as T 3 , we expect a 2% 
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error in the central temperature to give a 6% error in time. However, the stars that deplete lithium 
after (or while) contracting past the maximum temperature point always have t^est ^ *cont- Small 
corrections to the equation of state then have a large effect on the depletion time. Thus, our results 
are only reliable for stars that deplete lithium before reaching the maximum central temperature 
point, or that are younger than 



235 



M 



0.075Mt 







3 /2800ET\ 4 „ 

bar) Myr - 



(35) 



which limits us to masses M ^ O.O72M for lithium depletion by a factor of 100. The exact value of 
the low-mass cutoff is sensitive to the effective temperature, and Figure |3| displays the dependence 
of the cutoff mass on T e s and the level of depletion. For a given level of depletion, our calculations 
are valid only to the right of the corresponding curve in Figure ||. For masses lower than the above 
cutoff, one needs to include all of the relevant corrections to the equation of state in order to reliably 
predict lithium depletion at low levels and we refer the reader to the papers of CB97 and DM94 
for accurate calculations in this regime. 

In order to simplify the use of our results in evaluating observations and comparing to detailed 
calculations, we have constructed simple analytic fits to the solutions of equations (^J) and (|25|). 
The fits for age and central temperature of a star at a certain stage of 7 Li depletion (good for 
0.072 M £ M £ O.44M for 7 Li depletion by a factor of 100; the exact range is outlined in 
Figures |] and ||) are 
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(36) 
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3.06 x 10 G 



M 



0.149 



o.i. ",s / ^0.341 ( \ o.(«i 



x exp 



-0.0149 
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(37) 



K, 



where fj, = 4/(3 + 5X) for X = 0.7 and W = ln/ /ln/. The accuracy of these fits (as well as all 
the other fits presented below) is better than 1% for the age and 0.3% for the central temperature. 
The corresponding fits for 6 Li are 
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and 



T„ ( 6 Li) 



\O.lM y V3000AV \fioJ \]n2J v ; 



x exp 



•nil/ \5.7 / T \ 0.755 / ,, \ 9.55 / w \ 0.19 

-0.006159' °' 1M0 ^ ( TcS ^ ^ M ^ 



M J V3000Ky \nJ Vln2 



K. 



In the above equations the exponential factor accounts for the effects of degeneracy; for M > 0.2 Mq 
it can be set to 1 without sacrificing the accuracy of the fit. 

In our previous paper ( B 2 MU ), we compared our calculations for central temperature and age 
as functions of lithium depletion to those of CBP96. We found that the discrepancies between 
our results (using their effective temperatures) and theirs increased with the amount of lithium 
depletion. This difference has now been resolved. In the calculation of CBP96, the timestep did 
not scale with the rapidly decreasing lithium depletion timescale; as a result, when idest became 
much less than the Kelvin-Helmholtz timescale, the change in lithium abundance was no longer 
accurately computed. With the timestep adjusted to resolve idesti their code produces central 
temperatures (ages) consistent with ours to within 2% (10%) for 0.09 < (M/M Q ) < 0.2, and to 
within 3% (15%) for 0.2 < (M/M & ) < 0.35 (Baraffe 1997). Remaining discrepancies are now 
of comparable size to the differences between our results and those of DM94 (using the effective 
temperatures of DM94), and give us renewed confidence in our semi-analytic calculations. We also 
note that, since the time of elemental depletion is rather sensitive to the effective temperature 
(t oc 2^g 3 ' 5 ), the typical observational errors of ~ 5% in T e g lead to an uncertainty of ~ 15% in the 
age determination. On the other hand, the corresponding central temperature is rather insensitive 
to the errors in T cS flB^MUp . 



5.1.2. Comparison to Observations: Alpha Persei 
As an example of how our depletion calculations allow us to place constraints on ages of young 



clusters (first carried out in the Pleiades by Basri et al. 1996 and B^MU) we apply our work to 
the recently reported non-detection of lithium in faint objects in the a Per open cluster (Zapatero 
Osorio et al. 1996). We use their observed luminosities to infer the minimum age of the cluster, 
independent from the uncertainties in the effective temperature scale. From this minimum age and 
the observed effective temperatures we constrain the stellar masses. 



As described in B^MU, if lithium is undetected to some limit W = ln(/ //) in a star with 



—3.50 < log(L/L Q ) < —2.32 then the star must be older than 

/ i n -2.5 t \ °- 922 / ., x 2.52 / w \ 0.0769 
( _ M . 7 (a_V) p Myr . (40) 

This minimum age does not depend on the effective temperature; rather, the minimum arises be- 
cause electron-degeneracy pressure becomes important as low-mass stars contract flB^MU ). Equa- 



tion (EO) therefore describes an absolute minimum age for a lithium depleted star in the proper 
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luminosity range, independent of the uncertainties in the effective temperature scale. It is the 
dimmest star with no lithium that sets the most restrictive bound on the minimum age. Applying 
equation @ to AP 2790 (log(L/L ) = -2.49 ± 0.05), we find that the cluster, if coeval, must be 
older than 61 ± 7Myr, where the error comes from the uncertainty in the luminosity. The most 
recent age determination for this cluster from main-sequence turnoff by Meynet et al. (1993) gave 
53 Myr for moderate convective overshoot, and the conventional turnoff age without convective 
overshoot is 51 Myr ( flYlermilliod 1981] ). As with lithium dating of young stars in the Pleiades 
(Basri et al. 1996), we find that the nuclear age determination using lithium observations is con- 
sistent with traditional upper main-sequence fitting only if convective overshoot plays a significant 
role in the evolution of massive stars. The minimum age derived above implies that a Per is likely 
to be older than what has been previously inferred from main-sequence turnoff, even with moderate 
convective overshoot, in agreement with the conclusion of Zapatero Osorio et al. (1996). However, 
because of the uncertainty in the luminosity of AP 279, and if we allow a range of a few Myr in 
the ages of individual cluster stars, we cannot rule out the possibility that the lithium nuclear age 
and the conventional main-sequence turnoff age are consistent. If dimmer lithium-depleted stars 
are found in a Per, this consistency will be harder to maintain ( Basri 199 j ). 



Lithium ages of a Per and Pleiades are sensitive to the luminosities of their dimmest lithium- 
depleted members, tu oc L -0 - 9 (eq. |4Q| ), and, hence, to the assumed distance to the clusters. 



However, since the age determined from the main-sequence turnoff scales as £ms oc L ( Meynet 



et al. 1993 ), the difference between t\x and £ms is a weak function of the assumed distance. For 
example, the Pleiades would need to be at least a factor of 2 more distant for the lithium age to 
coincide with the upper main-sequence age with no convective overshoot. 

Given the inferred minimum age, we use the observed luminosity and lithium abundance to 
set a lower bound on each star's mass. Combining equations (||) and (^) yields the mass of a star 
that has a certain luminosity and age, 

" „. 122 f ' fV'f f«. (41) 



M Q ' VlO- 2 - 5 L y V 100 Myr 



J-cS 



Substituting the value of the minimum age into this equation gives the minimum mass for each 
cluster member that is fully convective and has not yet arrived on the main sequence. We illustrate 
this technique in Figure ^| for the stars AP279, AP272, and AP284. The shaded regions for each 
star represent the allowed values of mass and age given their level of lithium depletion. If we 
presume that they are all of an age equal to or larger than our minimum estimate, then we find the 
following minimum masses: AP143, O.3OM ; AP 296, O.25M ; AP284, O.18M ; AP268, O.11M ; 
AP 272, O.12M ; AP279, O.O9M ; and AP JO323+48530, O.O9M . It is important to keep in mind 



3 The latest spectral type a Per member considered by Zapatero Osorio et al. (1996) is AP J0323+4853. However, 
it is slightly brighter (log(L/ Lq) — —2.46 ± 0.06) than AP 279, and hence yields a less stringent minimum age. 



It has recently been reported (Martin & Zapatero Osorio 1997) that the rotation period of AP J0323+4853 is 



that the minimum mass determination, unlike the minimum age determination, depends on the 
effective temperature scale used, and so is subject to systematic errors. 



5.2. Beryllium Depletion 



Since our calculations require the star to be contracting and fully convective, or below the 
nuclear burning and radiative core lines in Figure [l], it is clear that our theory for 9 Be (standard 
rate) applies to stars with masses in the range 0.1-0.3 M & . All of these objects clearly deplete 
beryllium prior to reaching the main sequence, and for this mass range the age of the star at a 
given level of depletion is 
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Stars with M <^ 0.1M Q deplete beryllium while approaching the main sequence or residing on it. 
For these objects, we refer the reader to the results of Nelson et al. (1993) and CB97. For masses 
0.1M Q < M < 0.3M Q , our central temperatures (ages) differ by no more than 2% (13%) from 
those of CB97. Our comparisons are for abundances of more than 0.5 of the initial abundance, on 
account of the aforementioned error in the calculations of CB97. 

If the reaction 9 Be(p,a) 6 Li does proceed via a sub-threshold resonance (see § ||), then pa- 
rameterizing the cross section by the reduced width (see eq. Jl^j) implies that the age and central 
temperature at a given depletion level are 



100. 



O.IMp 



0.699 



3000 K\ 3 - 547 



x exp 



0.3435 



M ) V T cff 
O.1M \ 3 ' 612 / T cS 



Mo 



2.142 



M 



V3000K 



0.46 /,, \ 6.8 



W /0.1 ^ n 2 >7 

ry x 0.105 /q ^\ 0.217 

hi2 



(44) 
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7.6 hr. Using the formula for the radius (eq. ||), we estimate the breakup spin period (Ft w 2ty(R 3 /GM) 1/2 ) for the 
best guess at minimum age (61.2 Myr) and mass (O.O9M0) to be « 0.13 hr. Rotation is therefore unlikely to cause 
any deviations of the depletion time from our estimates (e.g., by decreasing the central temperature). 
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and 
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For example, a star of mass 0.2 M & at half-depletion will be 22% younger if the resonance exists 
and = 0.3. 



5.3. Boron Depletion 

Unlike lithium and beryllium, the higher Coulomb barrier of boron means that it is depleted 
when either hydrogen burning or radiative heat transport is important. In Figure ||, we show the 
central temperature at the time when one-half of 10 B (thin solid line) and n B (thin dotted line) 
are depleted, assuming that the stars are fully convective and not burning hydrogen. As one can 
see, this assumption holds only for a very narrow range of masses for 10 B, and does not hold for 
11 B at all. Thus we can only make qualitative conclusions about boron depletion during pre-main- 
sequence contraction. Only Nelson et al. (1993) and CB97 have tabulated the boron depletion. 
Neither of these authors displayed results for the less abundant isotope 10 B. If the resonance were 
not present in the p+ 10 B reaction (see § ||), then 10 B would deplete after n B and their boron results 
would have been in error. However, Figure [5| clearly shows that the resonance-corrected rate for 
this reaction is large enough so as to make it deplete before 11 B, in which case the aforementioned 
calculations are indeed valid. 

Our approximations are, however, applicable to stars in a certain mass range that deplete 
boron on the main sequence. While on the main sequence, the central temperature of a star is 
essentially fixed. Moreover, stars with M <^ 0.3M@ are fully convective on the main sequence. For 
such stars the formalism developed in § § is fully applicable, and, in particular, equations (p4|) and 
(|25| ) with a = const imply that the abundances of 10 B and 11 B will decay exponentially, at rates 
determined by the star's central density and temperature. Exponential decay of 11 B is consistent 
with its abundances at 1 and 10 Gyr in Table 2 of CB97; for stars of mass O.O85M £ M £ 0.13M Q , 
this depletion happens on the main sequence before the Hubble time. Moreover, the depletion time 
scale is an extremely strong function of mass within this mass range. 

It is clear that boron abundance in low-mass (M <^ 0.13M & ) main-sequence stars may provide 
a valuable constraint on their masses and ages, if boron can be observed. However, the strong 
dependence of depletion on the central temperature, as well as the high central density of these 
stars, necessitate sophisticated calculations like those of CB97 to accurately relate depletion time 
scales to stellar masses. 
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6. Conclusions 

The ability of stars to burn lithium, beryllium, and boron prior to reaching the main sequence 
has been recognized for quite some time. However, the observations of such depletion are just now 
becoming accessible for pre-main-sequence stars in nearby clusters. This has motivated our semi- 
analytic work on light element depletion for low-mass (M <J 0.5Mq) pre-main-sequence stars and 
the application of our results to observations of lithium-depleted stars. Lithium observations are 
proving to be an excellent age indicator in open clusters (see B^MU for a more extensive discussion) 
and our new minimium age determination for a Per (60 Myr) is in reasonable agreement with that 
found by upper main-sequence modeling with mild convective overshoot. However, if dimmer 
lithium-depleted stars are found in a Per, more overshoot may be necessary in order to reconcile 



the lithium nuclear age with the upper main-sequence turnoff one ( Basri 1998| ). In addition, the 
difference between the nuclear and upper main-sequence turnoff ages is rather insensitive to the 
distance to the cluster. We encourage observers to map out the lithium depletion region for pre- 
main-sequence stars in clusters of ages 10-200 Myr. This will provide an important check on the 
mass dependence of convective overshoot in the massive stars near the turnoff in these clusters. The 
main results of this paper are the fitting formulae for the age and central temperature as a function 
of mass, effective temperature, composition, and depletion level for 7 Li (eq. [Q and f37|), 6 Li (eq. 
|| and @), and 9 Be (eq. @ and @). These relations are robust and useful for evaluating 
observations and for comparing to detailed numerical calculations. Moreover, our results allow 
using the observed T e $ directly, instead of first placing a star on a theoretical evolutionary track, 
as well as exploring the effects of the different observational T e $ scales. 

Lithium dating of open clusters will most certainly help answer the question of how massive 
a star can become a white dwarf (WD). Ever since the white dwarf survey of young clusters 
of Romanishin & Angel (1980), there have been spectroscopic followups on the WD candidates 
(Reimers & Koester 1988). In addition to measuring the masses of the WD's, the ultimate goal 
of these efforts is to learn the range of stellar masses that can become isolated WD's. Crucial to 
this program is the ability to subtract the white dwarf cooling age from the cluster age so as to 
get the white dwarf birthday and thereby the mass of the progenitor star at that time. The most 
recent work of this nature is in NGC 2516 (Koester & Reimers 1993), which has an age infered 
from the upper main-sequence turnoff of 140 Myr (Meynet et al. 1993), giving an inferred mass of 
the progenitor star of around 6M Q (see also Reid 1996 for detailed discussion). Lithium dating 
of this (and other) young clusters would enhance our confidence in the extrapolation required to 
obtain the WD progenitor mass. 

In the cold atmosphere of a low-mass star, the typically observed absorption lines of beryllium, 
3130 and 3131 A, and boron, 2497 and 2498A flBoesgaard 1991 ), he in a part of the continuum where 



there are nearly no photons. Observations of pre-main-sequence depletion similar to those carried 
out for lithium thus await more sensitive instruments. There is, however, one intriguing possibility 
for diagnosing the elemental abundance of the surfaces of the low-mass main-sequence stars that 
might not have depleted beryllium or boron. Some of these low-mass objects are known to reside in 
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accreting binaries with white dwarfs, especially those with orbital periods below 2 hours (Warner 
1995). Most of these binaries are dwarf novae, in which the matter leaving the low-mass star due 
to Roche lobe overflow accumulates in the disk for some time before reaching a thermal instability 
strong enough to dump the matter onto the white dwarf. There have been recent HST observations 
of the white dwarf photosphere after these accretion events (Sion 1995; Sion et al. 1997; Cheng et 
al. 1997) which show that the WD is heated by accretion. Since the matter at the photosphere is 
actually from the low-mass star, there may be a chance of seeing some features from the beryllium 
and boron that are now residing in the WD photosphere, where the continuum in the blue/UV is 
much higher. This approach is complicated by the fact that nearly all binary evolution calculations 
suggest that low-mass objects in systems with such short periods have experienced substantial 
mass transfer. This means that the progenitors of these low-mass stars are likely to have been more 
massive in the past, and may have depleted beryllium and boron in their cores. 
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support as an Alfred P. Sloan Foundation fellow. 



-23- 



A. Fitting Formula for the Zero-Age Main Sequence 

Using the formalism presented in § we present a convenient fitting formula for the zero-age 
mainsequence (ZAMS) radius as a function of mass and effective temperature, 



R ni9of M Y 7m ( T ^ r 0-397 

= °- 129 [olMe J IsOOOK J 6Xp 



-0.192 



M \-5-73 / Tcff xl.48 



0.1M Q J V3000K 



(Al) 



This fit is good for 0.1 < M/M Q < 0.3 and 2500 K < T e{i < 3500 K. As a test of our approximations, 
we compare our radii to 1 Gyr isochrones from CB97 and DM94. We take the effective temperatures 
from the isochrone for each mass. Our radii deviate by no more than 2% (0.15 < M/Mq < 0.3) 
from Table 2 of CB97, and the maximum deviation for 0.1 < M/Mq < 0.3 is 10%. The increasing 
deviation at lower masses is due to the Coulomb correction to the equation of state. From Figure 
[l], the 1 Gyr isochrone of CB97 lies slightly above the maximum temperature reachable with our 
approximations (i.e., with Coulomb effects neglected). Our radii deviate no more than 8% (7%) 
from those inferred from Table 5 (6) of DM94. 

In Figure ^, we plot the surface gravity log g as a function of T e ^ for M/Mq = 0.12, 0.16, 
0.20, 0.24, and 0.30. Unfortunately, there is only one binary with stars in the mass range of 
our fit. This is the binary CM Draconis ( |Lacy 1977 ), which has a primary (secondary) of mass 



(0.237 ± O.O11)M [(0.207 ± 0.008)M Q ] and effective temperature of (3150 ± 100) K (same for the 
secondary). Using these parameters, we infer a radius of 0.246i? Q (O.221i? ) and a surface gravity 
log g = 5.03 (log g = 5.07); the values from Lacy (1977) are (0.252 ± 0.008)i? o [(0.235 ± 0.007)^©] 
and 5.01 ± 0.05 (same for the secondary). 

Similar fits were recently presented by Tout et al. (1996) for the ZAMS radii as functions of mass 
and metallicity. Our fits have the advantage of being in terms of mass and effective temperature, 
and so are independent of the model assumptions that are incorporated in mapping the metallicity 
to the effective temperature. 
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Table 1. Proton Capture Rate Parameters 



Reaction 


S 


a 


j 


References 


6 Li(p, 3 He) 4 He 


3.75 x 10 12 


84.13 


2/3 


1 


7 Li(p,a) 4 He 


7.20 x 10 10 


84.72 


2/3 


1,2 


9 Be + p (nonresonant)* 


2 x 2.18 x 10 13 


103.59 


2/3 


1 


9 Be(p, a) 6 Li (resonant) 


3.53 x W 16 9 2 


105.33 


1.367 


3 


10 B(p,a) 7 Be 


2.18 x 10 15 


120.62 


1/6 


1,2 


n B(p,a)2 4 He 


2.22 x 10 14 


120.95 


2/3 


1,4 



■"Includes contributions from 9 Be(p, 2 H)2 4 He and nonresonant 9 Be(p,a) 6 Li rates. 
References. — (1) CF88; (2) Raimann 1993; (3) Brown 1997; (4) Rauscher & Raimann 1996 
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Fig. 1. — Constraints on validity of Li, Li, and Be calculations. For Li we show the temperatures 
at which the abundance falls from one-half to one-hundredth its initial value (shaded region). For 
6 Li (dotted line) and 9 Be (dashed line) we plot only the curves along which one half of the element 
has been depleted. We also show the constraints on applicability of our depletion calculations. Stars 
above the line e rac j = e (heavy dot-dashed line) have radiative cores, while the curves L nnc = 0.05 L c 
and 0.50 L c (heavy solid lines) herald the arrival of the star on the main sequence. For comparison, 
we show (triangles) central temperatures of young (lGyr) main-sequence stars (CB97, Table 2). 
Coulomb corrections are unimportant except where the central temperature is near the maximum 
set by degeneracy (heavy dashed line), as described in § |5.1.1| , The maximum temperature curve is 
essentially a curve of constant degeneracy, as evidenced by the comparison with the curve where 
the electron pressure is 1.5 times its classical value (heavy dotted line). 
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Fig. 2. — Same as Figure [T], but for stellar age instead of central temperature. The main sequence, 
the maximum central temperature, and the constant degeneracy lines are not shown. The effective 
temperature scale T e ^{M) is from the results of Chabrier fc Baraffe 1997 . 
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Fig. 3. — Minimum mass, as a function of effective temperature, for which our 'Li depletion 
calculations are reliable. At a given effective temperature, we show the mass for which contraction 
at that T e g implies that lithium is depleted to the given amount ( 7 Li/ 7 Li Q = 0.01, 0.05, 0.15, 0.5, 
and 0.9) before the star reaches the maximum temperature set by degeneracy (eq. |]34|| ) . 
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Fig. 4. — An illustration of how the minimum age derived from equation ( f40| ) may be used in 
conjunction with the mass- age relationship, equation (|4l]), to set a lower bound on the masses of 
AP284, AP272, and AP279. We show the minimum cluster age (heavy vertical solid line) and its 
uncertainty (light vertical dotted lines). For each star, we plot two constant luminosity contours 
that correspond to the observational uncertainty in luminosity (light solid lines) and a contour of 
constant lithium abundance: AP284, [Li] = 12 + log(iV L i/iVH) < 0.23 (dot-long dashed line); AP272, 
[Li] < 1.94 (dashed line); and AP279, [Li] < 1.39 (dot-short dashed line). These constraints fence 
each star into a region of M-t space (shaded areas). The intersection of the shaded region with the 
minimum age curve sets a lower bound for the mass of each star. 
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Fig. 5. — The central temperatures at which the abundances of 10 B (light solid line) and n B 
(light dotted line) have decreased to one-half their initial values, presuming that the star is fully 
convective and the p-p energy generation is negligible. We also indicate (thick solid lines) where 
nuclear reactions contribute 5% and 50% of the stellar luminosity, and also (thick dot- dashed line) 
where a radiative region develops in the core. For comparison, we display (triangles) a 1 Gyr 
isochrone from Table 2 of Chabrier & Baraffe (1997). Clearly, any treatment of boron depletion 
must include hydrogen burning. For a small range of masses about O.1M , boron depletion occurs 
on the main sequence. 
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Fig. 6. — Each line is the logg-T e f[ relation for a low-mass star (of mass indicated on the plot) 
on the Zero- Age Main Sequence (ZAMS). We calculate the ZAMS using the method outlined in 
section 2.2. 



